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Inhomogeneous Patlak-Keller-Segel models and Aggregation 
Equations with Nonlinear Diffusion in M d 

Jacob Bedrossian* Nancy Rodriguez^ 



Abstract 



Aggregation equations and Patlak-Keller-Segel (PKS) models for chemotaxis with nonlinear 
diffusion are popular models for nonlocal aggregation phenomenon and are a source of a number 
of interesting mathematical problems in nonlinear PDE. The purpose of this work is twofold. 
C*~) . First, we continue our previous work [5], which focused on nonlocal aggregation, modeled with a 

convolution. The goal was to unify the local and global theory of these convolution-type models, 
p I i including the identification of a sharp critical mass; however, some cases involving unbounded 

| domains were left open. In particular, the biologically relevant case R 2 was not treated. In this 

paper, we present an alternative proof of local existence, which now applies to K d for all d > 2 
and give global results that were left open [5] . The proof departs from [5J [5] in that it uses a 
more direct and intuitive regularization that constructs approximate solutions on M. d instead of 
on sequences of bounded domains. Second, this work develops the local, subcritical, and small 
data critical theory for a variety of Patlak-Keller-Segel models with spatially varying diffusion 
and decay rate of the chcmo-attractant. 

> ' 

\q \ 1 Introduction 

In this paper we study several types of aggregation models with nonlinear diffusion and nonlocal 
self-attraction. The primary focus is to develop and extend the relevant local, subcritical, and 
small data critical/supercritical theory. These results exist in the perturbative regime as they all 
fundamentally treat the PDE as a nonlinear perturbation of the diffusion equation (see below for 
more information). Furthermore, we also present several non-perturbative global existence results 
for a class of critical problems as well. 

The first general class of systems we study are those where the nonlocal self-attraction arises as 
the result of a convolution operator 

u t + V • (uVK. * u) = AA(u), 

u(0,x) = u {x) e L 1 + (R d ; (l + \x\ 2 )dx)nL co (R d ), ^ ' 

where L+(M. d ;fi) : = {/ 6 L 1 (R d ;fi) : f > 0}. Equation ([T]) is of interest in mathematical biology 
as it models the competition between a species' desire to aggregate and to disperse. Dispersal is 
modeled via the, potentially nonlinear, diffusion term AA(u) and the aggregation is modeled via the 
nonlocal advective term V- (uVIC*u). The most well-known example of (H|) is the parabolic-elliptic 
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Patlak-Keller-Segel model, based on the original models of Patlak [37] and Keller and Segel [29] , 
For more information on the modeling aspects, see the [261 [25] for reviews of various chemotaxis 
models and [441 \TE\ l34"l [24l [To] for more general swarming and aggregation models. In this paper, 
we extend our recent work [5] to provide a more satisfactory and complete local and global the- 
ory for ([T]) on M. d for d > 2. In [5], we studied the local and global existence and uniqueness of 
bounded, integrable solutions to ([1]) in bounded domains for d > 2 and all space for d > 3. The 
primary goal of the previous work was to unify the existing Patlak-Keller-Segel global existence 
theory |41} |4"2| l4"3j [T4"l l3Tj \TT\ [TT] with the local existence and uniqueness theory for less singular 
versions of (pQ) [8]. In that work , the M? case was not treated due to several technical difficulties. 
Since this case is very important for mathematical biology, we make specific effort to treat this case 
and discuss the difficulties in more detail below. 



We present an alternative proof of local existence of ([I]) for a wide range of A and 1C which 
applies on M. d , d > 2, for solutions with bounded second moment. The new proof is based on a 
regularization directly on M. d , in contrast to [5]- One of the benefits of this regularization is that 
it allows one to rigorously justify the application of homogeneous Sobolev embeddings in formal 
arguments, which are crucial in deducing small data global existence and uniform boundedness 
in supercritical cases. We prove such results below following iteration techniques similar to those 
employed in [201 EJ E] . We also expand the global existence results of [5] to estimate the critical 
mass for kernels with a logarithmic singularity at the origin in M. d . 

The second class of nonlocal aggregation models we study is the variable-coefficient parabolic- 
elliptic Patlak-Keller-Segel model, 

u t + V • (uVc) = AA(u), 

-V ■ (o(x)Vc) + 7(z)c = u, (2) 
u(0,x) =u (x) e L 1 + (R d ;(l + \x\ 2 )dx)C]L co (R d ). 

We assume a(x) £ C 1 is strictly positive and j(x) £ L°° is non-negative. This PDE system is in 
many ways similar to (pQ), however they are not of the same form, as the solution to the equation 
for the chemical concentration, c, cannot be written in convolution form unless a(x) and j(x) are 
both constant. In this paper we develop the local, subcritical and small data critical/supercritical 
theory for ([2]). The proofs are analogous to those of ([T]) with additional complications arising due 
to the different nature of the estimates for c, Vc, D 2 c in terms of u (see Appendix). Analysis of the 
critical case and the identification of the critical mass has been completed by I. Kim and one of the 
authors in a separate work [3]. 

One of the most important properties of ([T]) and ([2]) is that each dissipate the following free 
energy 

T(u(t)) = S(u(t))-W(u(t)). (3) 
The entropy, S(u(t)), and the interaction energy (potential energy), W(u(t)), are given by 

S(u(t)) = / $(u(x,t))dx, 



W («(*)) = 9 / u(x,t)c(x,t)dx, 



2 



with c(x,t) = K, * u if the system is of convolution type. The entropy density (internal energy 
density), &(z), is a strictly convex function satisfying 




A' 



$'(1) = 0, $(0) = 0. 



(4) 



z 



In fact, both ([I]) and (|2|) are formally the gradient flows for ([3]) in the Euclidean Wasserstein metric 
(see e.g. PQ). For (pQ) and (J2j) there is no fully developed theory for making this precise; however, 
some aspects can be recovered and have proven very useful [10] - In any case, the free energy ([3|) 
plays an important role, especially in the global theory as in for example |4fi [T31 [T7\ [TT| \12\ 15]. 

Notation and Conventions 

We work on M. d for d > 2. For notational simplicity we denote the parabolic domain by = 
M. d x [0, T] and the standard L p norm by \\u\\ p = \\u\\ LP ^dy We also introduce the following notation 
for the k th moments: 



We let rj(x) G C£°(M. d ) with < r/(x) < 1 for \x\ < 1, r/(x) = for \x\ > 1, and r](x) = 1 for 
|x| < 1/2 be our canonical cut-off function and denote the standard mollifer J t v := rj * v. 

We use J\f to denote the Newtonian potential: 



Since we will be working with many largely irrelevant constants, we use the notation / < Pl k,... g to 
denote f < C(p, k, ..)<?, where C(p, k, ..) is a generic constant that depends on p, k etc. 

1.1 Definitions and Assumptions 

We consider the general class of kernels introduced in [S] , which includes fundamental solutions to 
elliptic PDEs and other commonly considered attractive kernels. 

Definition 1 (Admissible Kernel). We say a kernel K, G C 3 \ {0} is admissible if K, G W, ^ and the 
following holds: 

(KN) tC is radially symmetric, K,{x) = k{\x\) and A:(|x|) is non-increasing. 
(MN) k"{r) and k'(r)/r are monotone on r G (0,5) for some 5 > 0. 





By 'weighted Young's inequality' we mean for a, b > and 1 



i + i and e > 




(BD) \D 3 lC(x)\ < \x 



-d-l 
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This definition ensures that 1C is attractive, well-behaved at the origin, and has second deriva- 
tives that define bounded distributions on LP for 1 < p < oo. The obvious example of an admissible 
kernel is the Newtonian potential, which is effectively the most singular admissible kernel both at 
the origin and at infinity (in the sense that it decays the slowest). We remark that many of our 
results ( §2.11 and $3]) still hold if we replace condition (KN) with the assumption JC(x) = JC(—x), 
allowing for non-radially symmetric and general attractive/repulsive type kernels. 

We limit ourselves to diffusions that do not spread mass faster than linear diffusion; however, 
using the techniques of |17] one could also treat cases with fast diffusion. This is more general than 
the diffusion considered in [HI [5], which were restricted to degenerate diffusion. 

Definition 2 (Admissible Diffusion Functions). We say that the function A(u) is an admissible 
diffusion function if: 

(Dl) A G C^QO, oo)) with A'(z) > for z G (0, oo) and A(0) = 0. 

(D2) A'(z) > c for z > z c for some c, z c > 0. 

(D3) A'(z) < Ca when z < za for some za, Ca > 0. 

In particular (D3) implies A(z) = 0{z) as z — > and Jr u< m A(u) <c A ,R M for all R < oo. 

Following [U El [7] we use the following notion of weak solution, which is stronger than traditional 
distribution solutions. In d > 3, test functions are taken in H , whereas in d = 2 minor adjustments 
must be made, as discussed below and in [7J. By density arguments, this is basically the same as 
taking test functions in and requiring various regularity assumptions on the solution. However, 
but we prefer the current statement of the definition to emphasize the kind of test functions that 
we are interested in. Taking test functions in these spaces is important for the proof of uniqueness, 
which is based on an H" 1 stability estimate [El [2] . 

Definition 3 (Weak Solution in R d , d > 2). Let A and fC be admissible, and u G L^(R d ; (1 + 
\x\ 2 )dx) n L°°(R d ). If d > 3, a measurable function u : [0, T] x R d — > [0,oo) is a weak solution 
of flU) or © if u G L°°((0,T) x R d ) n L°°(0, T, L+(R d ; (1 + \x\ 2 )dx)), A{u) G L 2 {0,T,H 1 { 
nVc G L 2 (0,T;L 2 (M 2 )), u t G L 2 (0, T, H~ 1 (R d )), and for all test functions G L°°(0, T; H l { 

< u t , 4>{t) >fj-i X fji dt = - f f (VA(u) - u\7c) ■ V(j)(t)dxdt. 



H~ 1 xH 

JO J 

If d = 2, a measurable function u : [0, T] x R 2 — > [0, oo) is a weak solution of ([1]) or ([1]) if u G 
L oo ((0,T)xlR 2 )nL oo (0,r,L^(lR 2 ;(l + |2;| 2 )(ix)), VA(u) G L 2 (0, T, L 2 (R 2 )), nVc G L 2 (0, T; L 2 (R 2 )), 
u t G L 2 (0,T, V*(M 2 )), and for all test functions (f)(t) G L°°(0,T; V) we have, 

< u t , 4>(t) >v*xv dt = - [ [ (VA(u) - nVc) • V<p(t)dxdt, 
o Jo J 

where V = {/ G L°°(1R 2 ) : V/ G L 2 (R 2 )}. When solving CD) then c = JC * u and when solving 
then c(t) is the strong solution to —V • (a(x)Vc(t)) + j(x)c(t) = u(t) which vanishes at infinity. 

Remark 1. Due to the regularity we are imposing on our solution, we could equivalently require for 
all <f>e H 1 for a.e. t G (0, T) 



< Ut,(p>H-i X H^ 
with the modification in R 2 . 



J (VA(u) - nVc) • V<f>(t)dx, 
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Remark 2. The additional complication in R 2 is due to the fact that the norm ||V/||2 is not well- 
behaved in M 2 , since ||V/||2 scales like L°° in IR 2 . Indeed, there exists a sequence of Schwartz 
functions with ||V/ n ||2 = 1 and f n — > oo point-wise a.e. (consider f(x) = loglog(l + |a;| — 1 )l| a .| <;1 
and scaling f\(x) = f(Xx), A — > 0). In order to remove such pathologies from our space of test 
functions we follow [7], which also requires test functions to lie in L°°(IR 2 ). 

Following [5] , we now define a notion of criticality for (pQ) , which in general has no scaling sym- 
metries. However, a kind of scaling symmetry can be recovered in the limit of mass concentration, 
which in turn is expected to govern blow-up (see Theorem [3] below) . Suppose JC(x) ~ |a;| _c! / p as 
x — > 0, A(u) = u m with m > 1 and let u € n L\. Then if u\(x) = X~ d u(X~ 1 x) we have, 



lim J-(u\) = lim 
x^o A->0 



\d-dm r \~ d /P 

u m dx — / / u(x)u(y) \x — y\ /p dxdy 



m — 1 



From here we see the limit is independent of u unless d — dm = —d/p. If the limit is — oo we expect 
aggregation to dominate near mass concentration and if the limit is +oo we expect diffusion to 
dominate. As mass concentration should occur on vanishing length-scales, we may use this scaling 
heuristic to define a notion of criticality. The limiting case of p — > oo occurs when considering 
K,{x) ~ — log | X I clS X —7- and is discussed more below. Of course this corresponds to the 2D 
classical parabolic-elliptic PKS model, but we can consider the same singularity in higher dimensions 
(as done in for example [28]). 

Definition 4 (Critical Exponent). Suppose fC is admissible such that for some d/(d — 2) < p < oo 
(with the convention that p = oo if d = 2) we have D 2 fC(x) = O {\x\~ d ^ p ~ 2 ^ as x — > 0. Then the 
critical exponent associated to /C is given by 

v + 1 

1 < m * = £ < 2 - 2/d. 

p 

For the variable-coefficient Patlak-Keller-Segel system ([2]) we take m* = 2 — 2/d. 

Remark 3. Due to the monotonicity assumptions in Definition [1] (see also Lemma [1] above), for 
p < oo the definition is equivalent to requiring that KL{x) = 0(\x\~ d ^ p ) as x — > which is the same 
as requiring K, G L p ^ . Similarly, when m* = 1 we have that /C(x) = C(log \x\) as x — > 0. 

Remark 4. The variable coefficient system (J2]) should be roughly as singular as the constant coeffi- 
cient case, hence the corresponding definition in this case. 

Now we define the notion of criticality by relating the critical exponent of the kernel to the 
diffusion, again focusing on the limit of mass concentration. It is easier to define this notion in 
terms of the quantity A'(z), as opposed to using $(z) directly. This is not so surprising as A'(z) 
is precisely the local coefficient of diffusivity, and directly measures the strength of the diffusion 
relative to the mass density. 

Definition 5 (Criticality). We say that the problem is subcritical if 

limml — j-^r = oo, 

z— s-oo z m 

critical if 

A'(z) 

< liminf — )rzr < oo, 

z— too z m 
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and supercritical if 



Hminf^) =0 . 

z— >oo z m 



The following lemma, from [5], enumerates several important aspects of admissible kernels. Part 
(c) in particular provides a useful characterization of kernels with m* < 2 — 2/d. 

Lemma 1. Let K, be admissible. Then each of the following is true: 

(a) VIC e L^'°°, and if d > 3 then K, G L^'°°. Ifd = 2, then K, G BMO(R 2 ). 

(b) For all p, 1 < p < oo, 3 C{p) > such that for all f G L p (R d ), \\D 2 IC * f\\ p < C{p)\\f\\ p . 
Moreover, C(p) < p as p — > oo . 

(c,) Let 1 < m* < 2 - 2/d and /3 be such that 1 < P < d/2. Then, D 2 K, G Lf ' c °° j/ and on/?/ 
if VIC G Lf/ W - 1),0 °. 7/d > 3, tfien /C G L fWP-l),°° lf and only if D 2 K G L /3,oc flnd 
m* = 1 + 1/ p -2/d for some 1 < p < d/2 if and only if D 2 K G Lf ' 6 °°. JVote &y (Mj, £/ 
L> 2 /C G Lf o ' c °° /or some P > 1, then, D 2 )C G L^°° . 

We will also need the logarithmic Hardy-Littlewood-Sobolev inequality in order to relate the 
interaction energy to the Boltzmann entropy, as in for instance |214 [5] . 

Lemma 2 (Logarithmic Hardy-Littlewood-Sobolev inequality |18j). Let d > 2 and / G L\(R d ) be 
such that / log / G L 1 (M d ). Then, 



f(x)f(y)log\x-y\dxdy<^ f log fdx + C(||/||i). (5) 
Remark 5. One also has for all <5 > 0, 

log \x - y\ dxdy < Mi / /log/dx + (6) 

\x-y\<5 a JR d 

Also recall the well-known convolution inequality (see [32]): for all / G L p (R d ), g G L q (R d ) and 
K, G L''°° for 1 < p, q, t < oo satisfying 1/p + 1/q + 1/i = 2, 



f(x)g(y)K,(x - y)dxdy 



< P , q ,t,d ||/||p||flilg||£||L'>°°- (?) 



As noted above for d = 2 the admissible kernels are generally only in BMO(R 2 ) and hence can 
grow logarithmically at infinity. This introduces a number of complications for the local and global 
well-posedness. To begin with, in the proof of the energy dissipation inequality, one must ensure 
that the interaction energy of the approximate solutions converges to the interaction energy of the 
weak solution being constructed. However, JC * f will be unbounded for general / G L 1 n L°° hence 
more care must be taken than in d > 3. The dual of BMO is the Hardy space TL 1 @D], a strict 
subset of L 1 , which we define via duality, 

yi := sup 

KeBMO, || K || bmo = 



J )Cfdx. (8) 
Accordingly, we have the natural analogue of Holder's inequality |40| 



ICfdx 



<||/C||flM0||/||«i, (9) 
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which in particular implies K, * f G L°°(M. d ) whenever K, G BMO and / G rl 1 . The following lemma 
found in [2] provides sufficient conditions for / G L 1 such that / G % x and a useful estimate of the 
norm that supplies the convergence of the interaction energies. 

Lemma 3 ([2]). Let f G L l C\L p for some p > 1 and satisfy J fdx = 0, M\{f) < oo. T/ien f e H 1 
and 

WfWw <d, P \\f\\p + Mx(f). 

1.2 Statement of Results 

Theorem 1 (Local Existence and Energy Dissipation for Convolution- type Systems). Let d > 2, 
K, be admissble and uq G L\(M. d ; (1 + |x| 2 )cfa;) n L°°(IR d ). Then there exists a T > and a weak 
solution u(t) of Q which satisfies u(t) G C([0, T];L\(M. d ; (1 + \x\ 2 )dx)) n L°°((0,T) x M d )) and 
u(0) = iio- Moreover, J~(uq) < oo and it(f) satisfies the energy dissipation inequality, 

J c (u(t))+ [ f -^-\A'(u(s)Wu(s) -u(s)VJC*u(s)\ 2 dxds < JFIuq). (10) 
Jo J u { s ) 

For a rigorous interpretation of the free energy dissipation, 

/D[u]dx := f — I A'(u)X7u — uVIC * u\ 2 dx, 
J u 1 

see for example |19j . As mentioned above, the local existence for general models ([I]) for d > 3 was 
proved in [5] , which extended the existence results of |41|, [Til [8] . We present an alternative which 
is specialized to treating M. d that has certain advantages. In particular we treat M 2 . We also prove 
the corresponding theorem for variable-coefficient Patlak-Keller-Segel systems. 

Theorem 2 (Local Existence and Energy Dissipation for Variable-Coefficient Systems). Let d > 2, 
a(x) G C 1 be strictly positive such that a+ |Va| G L°° and let j(x) G L°° be non-negative. Ln d = 2, 
further suppose thatj(x) is strictly positive. If uq G L+(M. d ; (1 + |x| 2 )dx) n L co (U d ) then there exists 
a T > and a weak solution u(t) of ([2]) which satisfies u(t) G C([0, T]; L\(W d ; (1 + |x| 2 )cfa;)) n 
L°°((0,T) x M. )) and u(0) = uq. Moreover, J-(uq) < oo and u(t) satisfies the energy dissipation 
inequality 

TMt)) + [ [ — 7— r \A'(u(s))Vu(s) - u{s)Vc(s)\ 2 dxdt < FUq). (11) 
Jo J u{s) 

Uniqueness for convolution- type systems is proven in [5] for d > 3; the same proof works for 
([2]) using the elliptic estimates found in the Appendix. More recent work undertaken by J. Azzam 
and one of the authors [2] proves uniqueness for the d = 2 case. For completeness, we state a 
continuation theorem proved in [5] , which extends previous theorems stated in [271 [T7l [TT] . The 
extension to cover ([2]) is straightforward and is briefly discussed below in §3.11 

Theorem 3 (Continuation [5]). The weak solution to ^ or ^ has a maximal time interval of 
existence T+ = T^(uq) and if 

lim sup || (u — k)+\\ 2-m =0, (12) 
fc-»°°te[o,T*) 2 ~ m * 

thenT* = oo andu(t) G L°°((0, oo)xlR rf ). Herem is such that 1 <m< m* and\\va.'m.{ z ^ f00 A' {z)z l ~ m > 
0. In particular, if T* < oo, then for all p > (2 — m)/(2 — m*), 

lim HulL = oo. 



7 



One of the primary tools in the proofs of Theorems [I][3] is the use of Alikakos iteration methods 
commonly used in the study of these PDEs as for example \27\ [3"T| \T7\ [5j [3] . These methods are 
fundamentally perturbative in nature (as are the related methods of [381CL3]), depending on relatively 
crude Gagliardo-Nirenberg inequalities to overpower the nonlinear aggregation with diffusion only 
in certain regimes. In subcritical regimes this is sufficient and these methods prove global existence 
and uniform boundedness in L°°, as in |31|, I17j . In the critical and supercritical cases, one can prove 
the same provided that the initial condition is small in the corresponding critical norm and that 
the nonlinear diffusion compares favorably with the homogeneous diffusion even at low densities. 
As shown in for example [201 HH [381 E] , stronger decay estimates may also be deduced using various 
refinements of similar iteration methods. 

Theorem 4 (Subcritical and Small Data Theory). Let u G L\(W l ; (1 + \x\ 2 )dx) D (R d ) and let 
u(t) be the local-in-time weak solution to ([I]) or ([2]) with u(0) = uq. 

(i) (subcritical) Suppose that for z sufficiently large there exists some 5 > such that 

A'(z) > 5z m ~\ (13) 
for m > m* . Then the solution is global and u(t) G L°°((0, oo) x M. d ). 

(ii) (small data in critical and supercritical cases) Suppose instead (|13f) is satisfied for all z > 
for 1 < m < m* . If u(t) solves ([2]) and m < m* , then assume ■j(x) is strictly positive. Then 
there exists a constant eo = eo(8,m,d) > such that if 

\\UC\\\ 2-m < £(), 

2-m* 

then the solution is global (T±(uo) = oo) and u{t) G L°°((0, oo) x M. d ). Furthermore, if f{x) is 
not strictly positive and d>3 then there exists a constant cq = eo(<5, m, d) > such that if 

max ( ||«o||i, 1 1 no 1 1 2-m ) < eo, 

\ 2-m* / 

then the solution is global and u(t) G L°°((0, oo) x M. d ). 

Remark 6. Part (i) follows easily from the techniques used to prove the Theorem[3]found in |3 H ll7 [ lo] 
and will not be proved here. However, (ii) is not so immediate, especially for the variable-coefficient 
system ^ when j(x) is not strictly positive. 

Remark 7. One can check that if K. = M and A(u) = u m then if 1 < m < m*, ([T|) has a scaling 
symmetry which leaves the norm L^ 2 ~ m ^^ 2 ~ m *^ invariant. It is in this sense that the L^ 2 ~ m ^^ 2 ~ m *^ 
norm is critical. Hence, even when m < m*, part (ii) is still a small data critical result, which is 
expected. However the proof is more difficult in this case, as ([1]) has no controlled quantity which 
controls the L( 2 ~ m )/( 2_m *) a priori. 

Remark 8. We are unsure if the additional requirement on the mass in (ii) when j(x) is not strictly 
positive in supercritical problems is necessary. For the Alikakos iteration methods to work, the 
singularity of the advective term is quantified in terms of the elliptic estimates of Lemmas [11] and 
[T2l in the Appendix. In order to prove part (ii) for supercritical cases m < m*, one must use these 

2 — m 

estimates in a bootstrap argument to also control the critical norm L 2 ~ m * ; however, Lemma [TT] 
introduces a 'subcritical' lower-order term that seems to require the additional assumption in order 
to control. 
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For critical problems, deducing global bounds and decay estimates for larger data requires fully 
non-pertubative techniques that depend heavily on the energy dissipation inequalities (jlOp and 
(llip . The use of sharp functional inequalities to determine when mixed-sign energies such as ([3]) are 
coercive is the classical standard method of treating large data and determining the sharp critical 
mass, for example for the focusing nonlinear Schrodinger equations, marginal unstable thin film 
equations and, of course, PKS [35] HSJ [30], EJ H3 HU [5]. Our previous work of [5] treated the case 
m* > 1 for critical convolution- type systems in and in this work we complete the m* = 1 case, 
for example, now covering the variants of the classical Patlak-Keller-Segel system in M 2 . For Patlak- 
Keller-Segel, the proof that the critical mass is sharp follows easily from a standard Virial argument 
(see e.g. [14] [4T1 [43] \TT\ [5] or the more classical [35]), which can be modified in a straightforward 
manner to treat more general problems [5]. The corresponding program for the variable-coefficient 
systems ([2]) is a more difficult problem, completed by I. Kim and one of the authors elsewhere [4]. 
We also stress that (ii) is already known, e.g. [14\ 111] , but we restate and prove it to provide 
comparison with the other problems that do not have scaling symmetries. 

Theorem 5 (Global Existence in M. d , Critical case, m* = 1). Let A,]C be admissible and suppose 
uq £ L+(M. d ; (1 + |x| 2 )cfa;) n L°°(M. d ) and let u(t) be the local-in-time weak solution to ([I]) with 
u(0) = uq and define M = ||tto||i- 

(i) Suppose the problem is critical. Then there exists a critical mass M c such that if M < M c 
then u(t) exists globally. The estimate of M c is given below by (|14j) . If JC is bounded from 
below and 



I A'{z)z~ l dz < oo 

J 

then we may additionally assert u(t) € L°°((0, oo) x 



(ii) Suppose m* = 2 — 2/d, A(u) = u 2 2 l d and K. = N , the Newtonian potential. Furthermore, 
suppose that M < M c . Then u(t) exists globally and satisfies 



< 



(1+0 



Proposition 1. Let tC satisify IC(x) = —clog\x\ + o(log|x|), x — > for some c > 0. Then, the 
critical mass M c satisifies 

lim 4^ " ^-M c = 0. (14) 
z-^oo z log z 2d 

Remark 9. Note the additional requirements in part (i) of Theorem [5] in order to assert global 
boundedness. Although we believe that solutions are bounded (and in many cases likely decay as 
in (ii)) our proof cannot rule out an unbounded increase of entropy or interaction energy as the 
solution spreads without the additional assumptions. 



2 Local Existence and Energy Dissipation 
2.1 Local Existence for Convolution-type Systems 

This section focuses on the proof of local existence of weak solutions with bounded second moment. 
The proofs included here are simplifications of our work in [5] and combines techniques from the 
PKS model found in [14] and the non-singular aggregation-diffusion equations found in [8]. One 
advantage of the new techniques used here is the treatment of M 2 . Furthermore, in contrast to 



9 



the work in [5], we obtain solutions directly in R rf , removing the need for the intermediate step of 
finding solutions in bounded domains. Consider the regularized aggregation-diffusion equation 

u e t = AA e (u e )-V-(u € VIC e *u e ). (15) 

We define, 

A e (z) = [ a' e (z)dz, (16) 
J o 

where a' e (z) is a smooth function, such that A'(z) + e < a' t (z) < A'{z) + 2e. We denote, 

JCe (x) . = / I (S) M < 1 

1 K.(x) \x\ > 1 

Hence by Definition Q] /C e G C 3 (M d ). 

Proposition 2 (Local Existence for the Regularized System). Let e > be fixed and u${x) G 
L^(R d ; (1 + |x| 2 )cfe) n L°°(M d ). T/ien JT5J| fau a classical solution u e on [0,T] x M d /or all T > 

With U € (0) = J e UQ. 

To prove the above proposition we need some preliminary definitions. Define the Hilbert space 
V = {v G H 1 : \\v\\v < oo} , ||u||y = V< u j v >V> 

with the inner product defined via < u, v >v-=< u,v >#i + J |x| uvdx. Note by the Rellich- 
Khondrashov compactness theorem, V is compactly embedded in L 2 (M. d ) [22]. We will construct 
a weak solution to (|15[) with an analogous definition of weak solutions to Definition [3j For the 
remainder of the paper we denote the mass of the initial data, no, by M, i.e. ||«o||i = M. 

We prove Proposition [2] using the Schauder fixed point theorem, see e.g. [6]. The necessary 
compactness for the application is obtained via the Aubin-Lions Lemma |39| . Now we state and 
prove some a priori estimates that will be of used in the proof of Proposition [2j Some of these 
estimates are the same or closely related to estimates proved elsewhere (e.g. [HI [8l [5]). 

Lemma 4 (A priori bounds with linear advection). For fixed e > letu G L 2 (0, T; L 2 )nL°°(0, T; L 1 ) 
be given. Let u e be the global strong solution to 

u e t = AA e (u € )-V-(u € VIC e *u) (17) 
with initial data u = J e u (x) with u G L\(R d ; (1 + \x\ 2 )dx) fl L°°(R d ). Then, 

(i) ll u lL2(0,T;L2) < ^ 1/2 |klL-(0,T;L2) < T 1 / 2 \\ U 1| 2 exp (± || A/C e || 2 \\u\\ L 2 (p, T ,L*)T l/2 ) 

(ii) ||n e (t)||oo < ||uo||ooexp (HA/C'I^Hnll^^^)^/ 2 ). 

(in) ||V^ e (u e )||| 2(0)T;i . 2) < A^ll^llooJII^Hi + ||VX: e |||,oo||u||i«, ( o,r;ti)ll , * e lli~(o,r;£a)- 

(ivj ||Vu e || i a( 0i T,£a(R<J)) ^ VA e «) || L 2 (0iT ;i2) ■ 

(vj A^ 2 (^ e (t)) < C(M 2 (u ),e, ||«||L°°(o,T;i°°)) INU 2 (o,T;l, 2 )> IKI|i,r). 

(mj ||9t« e ||L2(0,T,H- 1 ) < C( e ; IM|loo( 0)T ..L°°)j \\u\\l^(0,T;L^), \\ u o\\l,T). 
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Proof. In what follows denote M := \\u e (t)\\i = ||«o||i- By (|17p and integration by parts, once on 
the diffusion term and twice on the aggregation term we obtain 



— J (u € ) 2 dx < -2 J |^(n e )| 2 + J V(u € ) 2 ■ V/C e * udx 

< - J V(u e ) 2 A)C e * udx 

< \\AJC e \\ 2 \\u\\2 J (u e ) 2 dx. 

Integrating implies, 

\\u(t)h < |K|| 2 exp Q||A/C e || 2 ||tt||i,i(o,T,L 2 )^ < Ikolb exp Q||A/C £ ||2||u|| L 2 (0iriL2) T 1/2 ^ . 

which gives (i). The bound (ii) follows similarly by estimating the growth of ||it(t)|| p and passing 
to the limit p — > oo. To continue, we need a bound on the L 1 norm of A e (u e ). Condition (D3) of 
Definition [2] implies that A e (z) < (Ca + 2e) z for some Ca > and sufficiently small z. Hence by 
Chebyshev's inequality, 

I A e (u)dx = [ A £ (u)dx+ [ A e (u)dx 

J J{u<l} J{u>l} 

< (C A + 2e)M + ^(||n|| 0O )A u (l) < (C A + 2e + A e (||n|| 00 ))M. 

We now turn to the less trivial estimate (Hi). Let tjji(x) := r](xR~ l ) for some R > 0, where 
rj is the smooth cut-off function defined above. Now take A = A e (u e )r]R as a test function in the 
definition of weak solution (Definition [3} , which implies, 

J (u t {t),Aj ds = J Q j (-VA e {u)+uVK, e *u)-VAdxds 



~~ Jo I VA ^ U ^ ' vA ^ dx + \ J Q T J\ VA( - U 



dxdt 



1 



T 



+ 2ll^ 7 ^ e lli oo ll^lli°°(0,T;L 1 ) y o J \ u<i \ 2 dx. 
Note also, we can apply the chain rule (see Lemma 14 in [5] or Lemma 6 in [8]) and get 

-J J(u € t (t), A)dxds < liminf J J A(s)dsdx. 

Furthermore, since A is increasing, we have J U A e (s)ds < A e (u)u. Therefore, 

V^ e (^)| 2 (o;)^ < ^(11^(0)11^)11^(0)11^! + IllV/C^lli^ll^lli^^^.^^ll^lll^o^.^^ 

+ [ T [ -VA e • V(A' m ) + \ |V(^)| 2 dx. 

JO JR/2<\x\<R Z 



- [ T [ \X7A e (v,Y 2 ' 
2 Jo J\x\<R/2 



>0 JR/2<\x\<R 

The latter error term can be bounded as follows 



f -VA*-V{A' m ) + \\V{A*T lR )\*dx< \ I \A\u')\ 2 \V m \ 2 dx 

JR/2<\x\<R Z Z JR/2<\x\<R 

<iT 2 ^(|K|U)PV)lli- 
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Therefore, since A € (u € ) G L 1 n L°°, we have for all T < oo, by taking R — > oo, 

l|V^ £ (u £ )||i 2(0jTiL2) < i(||n e || L -)||n e || Ll + ||V/C 6 |||o O ||5||| oo(0|T . £ i ) ||« 6 ||l. 

This completes the proof of (Hi). This bound, along with the fact that a' e > e, gives us the 
bound in (iv). 

The bound on the second moment of u e , (v), follows from the bound on the L 1 norm of A e (u): 
— j \x\ 2 u e dx = J \x\ 2 AA e {u e ) dx - j \x\ 2 V ■ (u £ V/C * u e ) dx 



\x\ 2 u e dx 



< 2 J x ■ VA e {u e ) dx + 2 J x- u e VK * u e dx 

<2d j A e (u e )dx + 2||V/C e * SllooIKH^ 2 (^j 

<2d j A e (u e )dx + 2M 1/2 \\VIC e \\ L ca\\u\\ L i (J \x\ 2 u e dx^j 
< 2 dJ A e (u e )dx + 2M 1 / 2 \\VIC e \\ L °o\\u\\ L1 (l + J \x\ 2 u € dx) 



1/2 
1/2 



An application of Gronwall's inequality then provides the quantitative result. Let <f> G C([0, T]; -ff 1 ), 
then by the definition of the weak solution (Definition [3l Remark [1]) and Cauchy-Schwarz, 

T pT 

{d t u"{t)A(t)) 2 dt< / ||VA(^)-^V/C e *5|| L2(Rd) ||V^)|| i2(Kd) ^. 
o Jo 

Therefore the bounds we have already obtained can be combined to imply \\dtu £ \\ L 2^ ^ T . H -i^ < 
C(T,u,u ). □ 

Now that we have all the required a priori estimates we are ready to prove Proposition [2j 

Proof. (Proposition [2]) Define the compact and convex subset of L 2 (0,T;L 2 ) 

S T = {v G L 2 (0,T;V) : |H| L2(0)T;H -i) + IM| l2 ( ,t ; v) < Co \\v\\l^(o,t-,l 1 ) < ||«o||i} , 

for some Co to be chosen below. Compactness in the L 2 (0,T; L 2 ) topology follows from the Lions- 
Aubin lemma (since V CC L 2 ) and the fact that St is closed due to the weak compactness of V 
and L 2 (0, T; H" 1 ), and the lower semi-continuity of the L 1 norm. 

Define the map J : St — > St as follows: given u G St we define Ju = u, where u is the solution 
to ifPT]) . Our goal is to apply the Schauder fixed point theorem. First, we verify that St is invariant 
under J. 



Step 1: (Invariance of St) Let u G St- Recall from clasical regularity theory (e.g. [33J or [22]) that 
since u satisfies the a priori estimate (ii) from Lemma [Hit is a global strong solution (of course it is 
not necessarily classical due to the potential irregularity of u). Furthermore, the additional bounds 
provided by Lemma H] plus conservation of mass give that u G St for C$ sufficiently large and T 
sufficiently small. Note that T and Cq depend only on M, the L°° norm of the initial data, and 
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e. Indeed, (i) provides a bound on |M|l2( 0j T;L 2 )i fa) provides a bound on ||Vti||£2( 0i T;L 2 )> an d (vi) 
provides a bound on \\dtu\\ L 2( ^ T . H -iy Moreover, the bound (ii) and (v) along with, 

1/2 

Ixl U 2 dx < WuWna ( / Ixf udx 1 1 1 1 1 ^ ^ 2 , 



provides an estimate on the first moment of the L 2 norm of u. Hence, J : St — > St- We are left to 
show that J is a continuous map in L 2 (0,T; L 2 ). 

Step 2: (Continuity in L 2 (0, T; L 2 )) We show J is continuous as a mapping from L 2 (0, T; L 2 ) 
to C([0, T\; H~ l ), which by interpolation against uniform bounds in H 1 provided by (iv), implies 
continuity in L 2 (0, T, L 2 (R d ) (since ll/IU < ||V/]|.!/ 2 ||/|| V 2 as can be seen from the Fourier trans- 
form). This approach, as opposed to working in L 2 directly, is convenient due to the nonlinear 
diffusion, which behaves most naturally in H . This argument is reminiscent of the H -1 stability 
estimate used to prove uniqueness of weak solutions [8j [7J [2] . Let {v n } n>0 C St be such that 
v n -> v in L 2 (0,T;L 2 (R d )). We show that Ju„ -»• in C([0,T];^- 1 (M d )). To this end, let 
n := — * ( Jw n — Jf ) and denote u n := Jw n and u := Ju. It is important to note that while the 
v n 's may not have constant L 1 norm, u Tl and u do since they have the same initial data. 

In order for us to proceed, we must show — = u n — u, in the sense of distributions, and 
that ||V0 n ||oo + llV^nlh ^c ,m,u ,t 1- To prove the former, it suffices to show that 4>(t) G L°°(]R d ), 
which implies — A(f> is a bounded distribution. All of these are proven by arguments found in the 
proof of Theorem 1 in [2] so we only briefly recall them here. In dimensions d > 3, these facts can 
all be established using the IP estimates on u n and u combined with Young's inequality and the 
fact that M e L d ^ d ~ 2 ^°° and V7V G L d l { - d ~ v >'°° . However, in d = 2 the problem is more delicate as 
M{x) = T^: log \x\ grows at infinity. Lemma [3J implies (p(t) G L°°, as f u n — udx = and u n ,u have 
uniformly bounded first moments. Proving ||V0(i)||2 ^ 1 can be shown using the Fourier transform: 
J u n — udx = implies u n — u(0, t) = and the uniform boundedness of the first moments of u n ,u 
implies Lipschitz continuity of u n — u. See [2] for more information. Now we may compute the 
evolution of HV^H 2 . = \\u n — u\\ 2 ^_ 1 : 

1 d f 

2jt / \ V ^)\ 2dX =< V<An(t),^V0„(t) >= - < <p(t),d t u n (t)-d t u(t) > . 

Therefore, using <f>{t) in the definition of weak solution, 

~ J \VMt fdx = - J (VA e (u n ) - VA e {u)) • V<j) n dx 

(u n — u)(V/C e * v) ■ V<j) n dx — j u(VIC e * (v n — v)) ■ V(f> n dx. 
h+h + h- 



/< 



We have dropped the time dependence for notational simplicity. Since A e is increasing, we have the 
desired monotonicity of the diffusion, 

{A € (u n ) - A e (u)) (u n - u)dx < 0. 



13 



Using integration by parts we have, 

A^ n V/C e * v ■ V(j) n dx 



1,3 1,3 

J 9 Xj 4) n d XjX .lC * vd Xi <p n dx - \®xi {(d X] 4>n) 2 ) d Xi IC € * vdx 
^ J d Xj (j) n d XjXi 1C * vd Xi (f> n dx + Y2 ^(dx^nfdx^lC 6 *vdx. 



t,3 hJ 



Hence, 



h < J \D 2 IC e *v\ \Vcp n \ 2 dx 
< ||D 2 /C e || 2 ||u|| 2 J \V<f) n \ 2 dx. 



Moreover, 



\h\ < ||w|| 2 ||VX; e * (w„ — w)||oo|| V^Ha 

^ IIV^IklK - ^||3/2||V0„|| 2 

S IK ~ HI^IK - u lll /3 ( 1 + l|V^n||l)- 

By the uniform bound on II^Hoo (by part (ii) above), the uniform bound on the mass in St and the 
regularization of /C, 



2^l|V0 n ||2 S \K - V\\r + (1 + ||u||i + \\v n - V f 2 ) ||V</> n || 2 . 

Integrating implies for some C > depending on the uniform bounds and e (using 0(0) = 0), 

||V0 n (i)||i < jf exp (cj\ + \\v(t')\\l + K(0 - v{t')\\ldA \\v n {s) - v(s)\\ 2 2 /3 ds. 

Since t < T we have, 

IIV^WIH < e OT +°Nli»(P.^,+cK-« ll^co, T .,, 2) / T K(a) _ , (s) ||f d , 

Jo 

By assumption, ||t>n(s) — v(s)||2 — > in L 2 ((0,T)) and hence also pointwise a.e. on (0, T). Since 

|| w n( s ) — ll^ 3 — 1 + ll u "( s ) — u (' s )ll2) by the dominated convergence theorem we have ||V0 n ||2 — * 
uniformly on (0, T). Therefore, on St, J is a continous mapping from L 2 (0, T; L 2 ) to C([0, T]; H -1 ) 
and hence by interpolation against the uniform H 1 bounds provided by (iv) in Lemma [H is also a 
continuous mapping from L 2 (0,T;L 2 ) to L 2 (0,T;L 2 ). 

Finally, we apply the Schauder fixed point theorem, which implies there exists a solution Ju = u 
with u € St for some T > 0. By the regularization of /C, it is straightforward to extend this solution 
to the regularized system indefinitely and to use a boot-strap argument to show that it is a classical 
solution to (fl~5l). □ 
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Proposition [2] provides a global family of classical solutions {u} €> q', however, we need to prove a 
priori bounds that are independent of the regularizing parameter e to prove obtain any compactness 
and pass to a solution of the original system (pQ) . The proof of compactness will strongly depend on 
the following iteration result. 

Lemma 5 (Iteration Lemma |31|, 117]). Let < T < oo and assume that there exists a c > and 
u c > such that A'(u) > c for all u > u c . Then if ||V/C * < C\ on [0,T] then \\u\\ < 
C2{C\) max{l, M, H^olloo} on the same time interval. 

The primary a priori bound required to obtain compactness is a bound on the L°° norm which 
is independent of e. The proof of Lemma 8 in [5] applies directly to the following lemma; however, 
we include the proof for the completeness. In this setting we simplify the proof by applying a 
homogeneous Gagliardo-Nirenberg inequality, Lemma fT3l 



Lemma 6 (L°° Bound of Solution). Let {u e } e>0 be the classical solutions to (|15f) on Dj-, with 
smooth, non-negative, and bounded initial data ^f e UQ. Then there exists C — C(||iio||i) H^olloo 

) and 

T = T(\\uo\\i, ||uo||p) f or an y P > d such that for all e > 0, 

IK(t)|| Loo(Rd) <C (18) 

for all t £ [0,T]. 

Proof. For notational simplicity we drop the e. The first step is to obtain an interval for which the 
LP norm of u is bounded for some p > d. Following the techniques commonly used, see for example 
[27J ED El (III E] > we define the function u e k = (u e — k) + , for k > 0. Due to conservation of mass the 
following inequality provides a bound for the LP norm of u given a bound on the LP norm of Uk, 



\uE < C(p)(k 



p-ii 



u\\i + \\Uk\ 



(19) 



We look at the time evolution of \\uk\ 



Step 1: 



< 



p J txg -1 V • (VA(u) - uVK 

—p(p — 1) J A'(u)Vufc • Vudx — p(p — 1) J uu^~ 2 Vtifc • V/C 
A'(u) 



4(p - 1) 



+ kp(p -1) J ul' 2 Vu k ■ V/C 
where we used the fact that for I > 

«(«*)' = (u k ) l+1 + ku l k . 
Integrating by parts once more and using Lemma [T] (b) gives, 



* u) dx 

* u dx. 

dx + p{p — 1) J u p k ~ L Vuk ■ V/C £ * u dx 
* u dx, 



(20) 



Jt \\u k \\ p p <-C(p) I A'{u) 



p/2 



dx — (p — 1) / u\MC e * udx — kp / u p k 1 A/C e * udx 



< -C(p)jA'(u) 

< -C(p)[A'(u) 



p/2 



Vu 



p/2 



dx + C(p) \\u k f p+l ||A/C £ * u\\ p+l + C(p)k \\u k \\ p p - L ||A/C e * 

: dx + c(p) (\\u k \\ p P X\ + ||u|iS5) + c( P )k (kh; + \\u\\i) . 
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Now, using (|19p we obtain that 

j t \\u h \\ p p dx<-C{p)jA\u) 



Vm 



p/2 



cte + C(p) + C(p, k) \\u k \\l + C(p, fc,M). 



An application of the Homogeneous Gagliardo-Nirenberg inequality gives that for any p such that 
d < 2(p + 1) (e.g. Lemma [T3] in the Appendix): 

a*! 



ip+i < 

lp+1 



a _> 



Vm 



p/2 



L- 



where ol\ = d/p, a 2 = (p+ 1) — d/2. From the the Young's inequality a r 6^ r ) < ra + (1 — r)6 (with 



5 || Vn p / 2 || i2 , Rd , and r = a±/2) we obtain 



ti 



ip+i < 1 



u 



+ r5 



lp+1 ~ "-up 

Above /3i, /?2 > 1. For /c large enough we have that A'{u) > c > over {u > k}; hence, if we choose 
5 small enough we obtain the final differential inequality: 

d 

<Yp° +U{p,k,rd) \\u k \ lf 



l < C(p,6) K||J 2 +C(p,k,r6) \\u k \\ p p + C(p,k, ||t*o[Ii). 



(21) 



The inequality (f2"Tj) in turns gives a T p = T(p) > such that is bounded on [0,T p ] and by 



()19[) then ||ti||p remains bounded on the same time interval. 

Step 2: We are not done since the bounds in (|2ip blow up as p — > oo. We follow the by now 
standard procedure and prove that the velocity field is bounded in L°°(M. d ) on some time interval 
[0, T]. Then this allows us to invoke Lemma [5] and obtain the desired bound. Since Vif G Lj and 
VK l R dv Bl (o) G for all q > d/(d — 1) (by LemmaCD (a)), then for all p > d/(d — 1) 

||V7C * n|| p < ||V/Cl Bl ( )||i||n|| p + ||V/Cl R d\ Bl ( ) \\ p M. 

By Lemma [T] (b) we also have, for all p, 1 < p < oo, 

\\D 2 K, * u\\ p < \\u\\ p . 

Then by Morrey's inequality we have V/C*n e G L°°(IR <: ') uniformly in e by choosing some p > d and 
invoking step one. Then Lemma [5] concludes the proof. □ 

Proof. (Theorem [1]) For all e > 0, let u e be the solution to [15] provided by Proposition [2j Lemma 
[6] provides a uniform-in-e a priori upper bound on ||it e (i)||oo on some time inverval (0, T) as well as 
a uniform bound on ||V/C * tt||oo- Hence, we may easily deduce the following a priori bound on the 
second moment: 

j t M 2 {u e (t)) = 2d J A e {u e )dx -2 J J u e x- VK e * u e dx 

<2d j A e {u e )dx + 2\\VK e *u e \\ OQ \\u\\ 1 1 /2 M 2 (u € (t)) 1 / 2 

< Ml J A^u')dx + (IKHoo + IMIO (1 + A^ 2 (n e (t))) . 

The first term is bounded in terms of ||u||i and || ti||oo as in the proof of Lemma HI By Gronwall's 
inequality we therefore have a uniform-in-e a priori upper bound on the second moment. Similarly, 
one may alter the proof of (Hi) in [4] to bound ||V-A||x^(o,T;L 2 ) independent of e using the uniform 
bound on 1 1 ^ e 1 1 oo ■ Using these a priori bounds we may follow the arguments of [5] (see also [H]) and 
prove the following lemma. 
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Lemma 7 (Local Pre-compactness in L 1 (My)). The sequence of solutions obtained via Proposition 
HI {u t } e> o, is pre-compact in L 1 ((0, T) x 



Equi-continuity relies on the fact that ||^4(ii e )|| i 2( 0i T;H 1 (Z))) < C uniformly in e. The proof de- 
pends on the domain size but when combined with tightness obtained from a priori bounds on 
the second moment (see below), this suffices. Hence, using Lemma [7] we may extract a subse- 
quence u tk — > u in L 1 ((0,T) x R d ). Additionally, since H'Wefclloo is uniformly bounded, by lower- 
semicontinuity, we also have u € L°°((0,T) x M. d ). Once the limit is extracted, we may further 
upgrade the convergence to C([0,T]; L 1 ) as in [HI [5] and hence to C([0,T]; L p ) for all 1 < p < oo 
due to the uniform L°° bound. The last remaining technical point is to ensure that the limit is 
indeed a weak solution of the original system in the sense of Definition O This may be done using 
a limiting argument and the a priori bounds available on u as in [5j [8] . 

Proving the energy dissipation inequality in M? also presents an additional complication. As K, 
is potentially logarithmically unbounded at infinity and hence the convergence of the interaction 
energy is non-trivial, we cannot follow the exact argument from [5]. Hence, let u(t) be the unique 
weak solution constructed above and {u e } e>0 be the solutions to the regularized system (fT5|) . We 
may follow the proof in [5] to prove convergence of the entropy and the energy dissipation, hence, 
we need only focus on the interaction energy. To this end, 



u e (t))C e * u e (t)dx - J u{t)K * u(t)dx = I u(JC e -IC)*u + J ulC * (u e - u)dx 

+ J(v e - u)K, e * u e dx 
■ - Tl + T2 + T3. 

Since JC e (x) = tC(x) for all \x\ > 1 we have, 

|T1| < UnldlK/C 6 -K) *u\\ OD 

< IMIilMloo||£ e - m^iB^o)) o. 

Now consider T2. By the duality of BMO and T-L 1 and Lemma [3] we have, 

\T2\ < ||u||i||/C e * (u € -u)||oo 

< \\IC e \\BMo\\u e ~ uWu 1 

< \\u e — u\\ p + J \x\ \u € — u\ dx, 

for any 1 < p < oo. However, u e — > u in C([0, T]; L p ) for all such p, so the first term is not an issue. 
To deal with the second term we first use Cauchy-Schwarz 

1/2 , . \ 1/2 

\x\ \u e — u\ dx < I / \x\ \u e — u\ dx } / \u e — u\ dx 



However, since u € — > u in C([0, T]; L 1 ) and both u e and u have uniformly bounded second moments 
on [0, T] we have that T2 — > 0. The final term, T3, follows similarly. Hence, the energy dissipation 
inequality holds in ]R 2 . □ 
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2.2 Local Existence for Spatially Inhomogeneous Patlak-Keller-Segel Systems 

In this section we prove local existence of solutions of ([2]). Once again, we begin with a regulariza- 
tion, which in this case is given by 

u% + V • {u e VJ e c e ) = AA e (u e ), 

-V • (a(x)Vc E ) + -y(x)c e = J e u e , (22) 
u e (0,x) =J e u Q (x), n G L 1 + (R d ;(l + \x\ 2 )dx)nL°°(R d ), 

where A e (u £ ) is the parabolic regularization given by (|16p . This particular regularization is conve- 
nient because (I22j) satisfies a natural energy dissipation inequality, 



F(u e ) = j § e (u e )dx -^J u e J e c e dx. 



Indeed, using 



one can formally obtain 



ld_ 

2dt 



J u e J t c dx = j u e t J e c dx, 



^-F e {u e ) = - I — \(A e y(u e )Vu e -u e VJ e c\ 2 dx. 
dt J u e 



Moreover, the regularization of the chemo-attractant only occurs on the right-hand side, which 
proves to be of some technical use. 

The proof of Theorem [2] is very similar to that of Theorem [TJ one simply needs to check that c 
satisfies properties sufficiently similar to those satisfied by convolution-type systems. Most of these 
properties are obatained via standard elliptic estimates, which we state and prove in the Appendix 
to be referenced as needed. As was done for the convolution-type systems, we first need equivalent 
estimates as those provided by Lemma [4] for the system 

4 + v • (u e VJ e c) = AA e (u e ), 

-V • (a(x)Vc £ ) + j(x)c e = J e u, { ' 

where u € L 2 (0, T; L 2 ) n L°°(0, T; L 1 ) is given. As above, this will allow us to obtain the necessary 
compactness and a priori estimates to apply the Schauder fixed point theorem and obtain a family 
of classical solutions for (|23p , as in Proposition [2j 

Proof. (Theorem [2]) As this proof follows that of Theorem (pQ) we simply mention and work out 
the differences. The first step is to prove the existence of solutions to the regularized system (I23p . 
This is done, as in the proof of Proposition ([2]), by a fixed point argument. Hence, we need to show 
that for a given u the system ()23p has a solutions that satisfies the appropriate bounds. Examining 
the proof of Lemma ([3]) we see that all the bounds hold provided we have the bounds on || A^cHoo 
and llV^clloo. Using Young's inequality for convolutions, when 7 > we obtain 

IIAJ^cHoo < e ||c|| 2 < 7 |H|2 
IIVJ^clloo < e ||c|| 2 < 7 ||£t|| 2 . 

The estimates are different in the case when 7 > 0, as we only have estimate (|36p (necessarily 
we are assuming d > 3). We instead have, 

HAJ^Hoo < e ||c||_M_ < ||u|| 6 /5 < INl2 /3 |Nll /3 - 
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Although this estimate involves a lower power of ||u||2j this does not pose a problem for our needs. 
A similar bound for ||Vi7 £ c||oo can be obtained. These bounds are sufficient to obtain a family of 
solutions on {u € } e>0 on [0, T € ) by following the proof of Proposition [2j Next, we need to prove a 
uniform- in-e L°° bound of the solutions, similar to Lemma [6] Recall that the first step in the proof 
of Lemma [6] was to obtain a bound on ||tt e || p for p > 1. This required a bound of the form 

||Ac e || p < C(p,a,d)\\u e \\ p , (24) 

for example, to obtain the inequalities after (|20]h If 7(x) is strictly positive, a proof of (|24p can be 
found in Lemma [12] in the Appendix. For the remaining cases, d > 3, this estimate is replaced by 

||Ac e ||p < C(p,a,d) (\\u e \\ p + K||i), 

which is obtained from Lemma [TTJ noting that 2 p+d ^ P anc ^ interpolating between the LP and 
L 1 norms. The the next step was to use the bounds on the LP norms of u e to deduce a bound 
on ||Vc e ||oo, independent of e. By Morrey's inequality and then the Gagliardo-Nirenberg inequality 

ll Vc l|oo < ||^> 2 c||p + ||Vc|| p 



< \\D 2 c\L + 



c 



If 7 is strictly positive then we may use Lemma [12] and we have 

llVdloo < 



If 7 is not strictly positive we have to use the LP estimate (|36p and Lemma [IT] which holds in d > 3 
for p > d/(d - 2): 

||Vc||oo < a ,p,d \\u\\ p + \\u\\_pd_ < \\u\\ p + ||n||i. 

d+2 

Hence, we may proceed as in Lemma [6] and apply Lemma [5] and deduce the requisite local in time, 
uniform-in-e L°° estimates. 

The last primary estimate we need before we can obtain enough pre-compactness to extract a 
subsequence which converges towards a solution to the original model ([2]) is a uniform estimate on 
||j4(m)||£2(o i t,.H' 1 )- Examining the proof of estimate {Hi) in Lemma H] one can see that the bound 
on HVclloo is sufficient. Analogously as to above, pre-compactness of {u e } e>0 follows from these 
uniform estimates. Finally, to obtain a solution to ([2]) the only real difference with convolution-type 
systems is that we need to check that ||it e Vc e — itVc||x,2 — > as e — > 0. Indeed, we know that 

|| u e Vc e — uVc\\l2 < KHl 00 II V(c e — c) 1 1^,2 + HVclloolItt 6 — u\\ L 2. 

Furthermore, again using (I38p we have 

||V(c e -c)|| L 2 <||£>V-c)||l2 + || C £ - C || l2 . 

The result follows from Lemma [12] or Lemma [TT1 and the convergence in LP of J e vf to u. Hence, 
one may deduce that u is a weak solution to ([2]) in the sense of Definition [3J We are left to prove 
the energy dissipation inequality (jlip . Convergence of the entropy term follows exactly as the proof 
of inequality (|10p . given in Theorem [T] The interaction energy on the other hand requires a little 
more. As above, 

\\J e U € C e - UC\\ L 1 < ||J" e 1i e ||l||c e - c|| L oo + ll^lloollJ,^ - U||jr,l. 
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Given that u e — > u in L 1 and that we have uniform bounds on llc 6 ^ by and = M we only 

have to verify the uniform convergence of c e to c. It follows from similar arguments as above, using 
Morrey's inequality and (13811 

l|c e - c|| L oc < ||V(c e - c)\\ p + \\c e - c\\ p 
<\\D 2 (c £ -c)\\ p + \\c e -c\\ p , 

Hence, the result follows from Lemma[TT]and ()36|) or Lemma[12]and the LP uniform-in-time conver- 
gence of u e to u. The last step that requires attention is proving the convergence of the free energy 
dissipation (also called the generalized Fisher information), 



D[u(s)]ds = [ [ —, r \A'(u(x, s))Vu(x, s) — u(x, s)Vc(x, s)\ 2 dxds. 

u[x, s) 1 



However, as in the convolution-type case in [3 [5], it follows from a weak lower-semicontinuity result 
due to Otto [36] (see also US]). □ 



3 Continuation, Subcritical and Small Data Theory 

3.1 Continuation for Spatially Inhomogeneous Patlak-Keller-Segel Systems 

The proof of Theorem [3] is very similar to that of Lemma [6] and largely follows the work of \17\ [TT] . 
Both of these results quantify the strength of the nonlocal advection term with the inequality 

l|A/C * u\\ p < Ptd \\u\\ p . 

For the system ([2]), unless j(x) is strictly positive, we only have the estimate provided by Lemma 
[TTlin the Appendix: for p > d/(d — 2), 

ll Ac llp ^p-d \\u\\p + \\u\\ pd < \\u\\ p + \\u\\x. 

2p+d 

The lower-order term (which is constant due to conservation of mass) can be safely added into 
other lower-order terms in the argument which are already present due to the inequality (|19p . The 
difference between Lemma [6] and Theorem [3] is the use of a better Gagliardo-Nirenberg inequality 
to control the highest order contribution from the advection term: 

\\u k \\ p+ i ^ Ffcll^^^m ) J vn fc ax, 

with a.2 > determined via scale invariance and homogeneity (Lemma I13p. The latter factor can 
be related to the nonlinear diffusion and the former factor is the quantity appearing in <\12h . The 
key point is that the condition (|12j) then implies we may make the highest order contribution of the 
nonlinear advection arbitrarily small compared to the diffusion by choosing k large. 



3.2 Subcritical and Small Data Theory 

Proof. (Theorem [4]) We proceed with the formal computations, noting that the computations are 
completely rigorous by appealing to the parabolic regularization above. The proof follows a similar 
outline to the proof of Lemma [6] or Theorem [3l which is now standard in the literature on PKS 
and related models. As remarked above, the subcritical case is a standard variant of the proof of 
Theorem [3l |31 | \T7 \ [5] : therefore, we only prove (ii). We first prove the result for convolution- type 
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systems and for variable-coefficient systems where 7 is strictly positive. 



Let q = max((2— m)/(2— m*), 1). We have two cases to consider: m* = 2 — 2/d and m* < 2—2/d. 
We make this distinction in order to take advantage of the additional regularity provided by (c) of 
Lemma [T] when m* < 2 — 2/d; hence, in this case we let j3 be such that m* = 1 + 1//3 — 2/d for 
some 1 < (3 < d/2. Note that by Lemma \T\D 2 K, G Z/>°°. 



Step 1 (LP bound): Choose p such that if m* = 2 — 2/d then q < p < 00 else q < p < {3/({3 — 1) 
Computing the time evolution of ||it(i)|| p making use of ()13j) we obtain, 



j t \\u\\* < -4(p - 1)5 



m—l 



VvP' 2 



dx — (p — 1) / u p Acdx 



If m* = 2 — 2/d, we may use Holder's inequality and then Lemma [JJ for convolution- type systems 
or Lemma [12] for variable coefficient systems to obtain, 



vPAcdx 



< P ,tc \\u\\p+v 



On the other hand, if the system is of convolution-type and /3 > 1, since D 2 1C £ LP* 00 we can apply 



u p AK, * udx 



^ ll U llQ P ll U l|t|l A ^llL«9,->, 



with the scaling condition 1/a + 1/t + l/j3 = 2. Choosing t = ap implies, 



1 



a 



l/P 



l + l/p 



(25) 



Notice that from our choice of p, 1 < 1/p + 1//?; thus, 1/a < 1. When m* = 2 — 2/d, we define 
t = ap = p + 1. By the Gagliardo-Nirenberg inequality Lemma [TBI using m < m* we have, 



with 



and 



Map ^ \\u U q 



2d 

a\ = — 



Vu (p+m_1)/2 
(p-q/a) 



\ on/2 

dx 



(26) 



p \q{2 — d) + dp + d(m — 1) J ' 



a 2 = 1 - «i(p + m - l)/2. 
Notice that by our choices, a\{p + l)/2 = 1. Hence ([26]) implies, 



|ii-n?<(^)ikiir {p+1) -^) 5 ) 



.m— 1 



v« p / 2 



(27) 



€ L°°(IR+). If m < 
is not conserved. However, in that case (|27p also holds for p = q, and 



If m = m* then q = 1 and for mass sufficiently small we have \\u 
then g > 1 and ||u() 

therefore if ||^o||g is sufficiently small, (I27j) implies ||u(t)||q < ||uo||g for all t > 0. Hence it follows 
again from (|27p that there exists a constant C(p) such that if H^oHg < C(p) then ||u(t)|| p G L°°(1R + ). 
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Step 2 (L°° bound): Once more we use Lemma[5j which implies that if Vc 6 L°°((0, oo) x then 
u e L°°((0,oo) x If m* = 2- 2/d, by Lemma □ or Lemma [El 

ll^ 2 c|U+i < IMU+i- 

For systems of convolution type we have by weak Young's inequality, 

||V/C * 1i || 2d/ (<2-l) < ||V/C|| £ ,d/(<J-i),<x.||«||(2d-l)/(d-l) ) 

and for variable coefficient systems we may use (138p and Lemma [T2l 

||Vc||d+i < ||I? 2 c|| d+ i + ||c|| d+ i < ||u||d+i. 

By Step 1, we may choose H^oll? sufficiently small such that both of these norms are uniformly 
bounded, which is sufficient to apply Morrey's inequality and conclude that Vc 6 L°°((0, oo) x 
In the case of convolution- type systems with 1 < m* < 2 — 2/d we may only bound L p norms of u 
uniformly in time for 1 < p < j3/((3— 1). Since V/C £ Lj oc and V/Cl^ye^o) G L q for all q > d/ (d— 1) 
(by Lemma [1]), we have for any q > d/{d — 1) 



I V/C * u|| g < ||V/Cl Bl ( )||i||n|| g + ||V/Clj 



\Bi(0)ll9 



M. 



Since m* < 2 — 2/d, (3 > 1, and then we may choose g € (<i/ (<i — — 1)], since in this case 

necessarily d > 3. By the weak Young's inequality we have, 



\D 2 K. * u\\d+i < ||£> 2 /C|| L e,oo \\u\\ /3(d+i) , 

d(/3-l)+2,3-l 



where 



1 < 



< 



/8 



d(/3 - 1) + 2/3 - 1 ~ /3 - 1 ' 



and ||«| 



By Step 1, we may again choose ||ito||n sufficiently small to ensure that \\u\ 

d([i-l)+2f}-\ 

are both uniformly bounded. Therefore, we may again apply Morrey's inequality and deduce 
V/C * u £ L°°((0, oo x R d ) and conclude the proof using Lemma [SJ 

In order to treat the case d > 3 and 7(2;) not strictly positive we use the less optimal estimate 
provided by Lemma [TT1 Hence beginning as before we have for p > q and p > 1, 



d 



dt"~ 

Using Lemma [TT] we have, 



u\\p < -4(p - 1)5 y u™" 1 Vu p / 2 2 dx - (p - 1) J u p Acdx. 



u p Acdx 



11 



ip+i 



p+i + IrIIp+i ll^ll Cg+jjj ' 



2p+2+d 



(28) 



note that necessarily p + 1 > 2 > d/(d — 2) and 2p+2+rf — P + 1- By the Gagliardo-Nirenberg 
inequality, Lemma [T5| 

2 \ a i/ 2 



with 



Mlp+l ~ IKlg 



2rf 



Q2 



ai = — - 



V-u (p+m_1)/2 
(p - pq/(p + 1)) 



dx 



(29) 



p \g(2 - d) + dp + d(m - 1) 
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and 

02 = 1 — ai(p + m — l)/2. 
As above + l)/2 = 1. Let 6 £ (0, 1) such that \\u\\ (p+ i )d < Af 1_e ||u||p +1 . Therefore by fl28|) we 



2p+2+ 

have, 



If q = 1 we may choose M sufficiently small (depending on p) so that the first term is negative (in 
fact, as negative as we require, say < —5/2). Since again by interpolation \\u\\p < + M we 

then have 

^u\\p<- 5 -\\u\\p + C(p)M^ 1 + 6 -M. 

Hence for all p we may choose M sufficiently small such that the LP norm of u(t) is uniformly 
bounded in time. From here we may proceed similarly to above (using Lemma [11] instead of Lemma 
n~2|) and eventually invoke Lemma [5j If q > 1 more thought is necessary since ||it(t)||g is not a 
conserved quantity and the continuity argument used for convolution-type systems will need to be 
refined due to the presence of the lower-order term. From the previous argument in the q = 1 case 
it is clear that the result will follow if for any e > 0, we may choose ||wo||$ sufficiently small to ensure 
||g < e uniformly. To this end we note that (1301) holds for all q: 

The result will follow from a continuity argument. If ||uo||q is chosen such that 

II ""llg 



then for some time (depending on M) this implies, 

d „ „f, . 5 , 



which in turn implies, 



, t H\l<- U -\\u\\lll + C( P )M^ 

<- 5 -\\u\\l + C(p)M^ l + 6 -M, 



\u(t)\\l < max (\\u\\l I (c{p)M« +l + 5 -m\\ 



From here it is clear that if the mass is chosen sufficiently small then ||it(t)||q can be bounded 
uniformly by a sufficiently small constant. Hence we may proceed as after (I30p above to finish the 
proof. □ 



23 



4 Large Data Theory for Critical Convolution-type Systems 

As in [5] 1141 ITT], global existence in critical cases will follow from continuation (Theorem [3|) in [5] 
provided Js u>k \ ulogudx is bounded for large k, which will be deduced using the energy dissipation 
inequality (llOj) . There are two difficulties here. Firstly, due to (D3), the nonlinear entropy J &(u)dx 
and the Boltzmann entropy J u log udx are no longer uniformly bounded from below as the solution 
spreads out. Moreover, and perhaps more fundamentally, is the possibility of /C not being bounded 
from below, which again can cause the potential energy to grow unboundedly as the solution spreads 
out. Both difficulties will be overcome by the following two standard lemmas. The first shows that 
the decay of the entropy is bounded from below by controllable quantities. 

Lemma 8 (Entropy Lower Bound). Let A be admissible, u £ L+(M. d ; (1 + |x| 2 )<ix) n L°°(R. d ) with 
M = ||u||i. Then for all e > 0, 

(e d l 2 M\ 

ulogudx > Mlog [ _ d/2 j - eM 2 (u), (31) 



7T 



and 

e d l 2 M\ 



J <$>{u)dx > C A M log — ^- - eC A M 2 (u). (32) 



Proof. Following [13] . by Jensen's inequality for probability measures (dfi = jfadx), 



I 



u(x) log u(x)dx + e J \x\ 2 u{x)dx = — M J log I — ^— — j j^dx 



[ f e~^dx 



M log 



' e d ' 2 M \ 

n d/2 J 



Therefore, © holds. Next, by (D3), for some 6 > 0, A'(z) < C A for z < 5. Let h(z) 
Jj 2 A' (s)s~ 1 ds and note that f &(u)dx = f Rd Jq h(z)dzdx. For z < 1 we have, 

h ( z ) = - f 1 ^l d s 



s 

j 1 ^l ds _ Ca t -ds 

min(z,<5) ^ 



S 

> -C + C A logz. 

Therefore, since logz is integrable at zero, we have the following by Chebyshev's inequality, 
J &(u)dx = J J h(z)dzdx 

/* f'U f'l 

> / l{ u <i} / h(z)dz + l{ u >i} / h(z)dzdx 
J Jo Jo 

> J 1 {u<i} (C A (ulogu - u) - Cu) - Cl {u>1} dx 
>C A [ u log udx-C(M). 

J{u<l} 

Therefore, using ([3T|) we obtain ([32]) . □ 
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The following lemma establishes a uniform bound on the second moments for critical problems 
with to* = 1. 

Lemma 9 (Second Moment Estimate). Let A and K, be admissible and critical with to* = 1. Then, 

M 2 {t) < M 2 (0) + M (Ci + C 2 M) t 

for some constants Ci > 0. 

Proof. We argue formally, noting that the computations can easily be made rigorous with standard 
arguments. Computing the time evolution of the second moment, 



d 
di 



M 2 (t) = 2d j A{u)dx + j j u(x)u(y)(x - y) ■ V/C(x - y)dxdy. 



By (Dl) and (D3) of Definition [21 and criticality in the sense of Definition [5l we necessarily have 
A'(z) < 1, and hence J A(u)dx < M. By Definition Q] and m* = 1 we have, 

\(x-y)-VK(x-y)\ < C. 

Therefore by integration, the lemma follows. □ 

We now prove Theorem [5l 

Proof. (Theorem [5]) We proceed formally, noting that the arguments can be made rigorous with 
the regularization procedure of the local existence theory. We begin by proving (i). We first prove 
the result for the case when either K is not bounded below or the diffusion is not degenerate. As 
noted above, we will not get uniform-in-time bounds. Recall the energy dissipation inequality (jlOp . 



J §{u)dx - i j ulC* udx < F(uq) := F . 



By the asymptotic expansion of the kernel assumed in Proposition [T] and (BD), we have that for 
all e > 0, 3 5, R > such that, 

§(u)dx H ^— / / u{x)u(y) log \x — y\ dxdy 

2 J J\x-y\<8 

< Fq H — / / u(x)u{y)K,{\x — y\)dxdy — C j j u(x)u(y) log \x — y\ dxdy. 

2 J J5<\x-y\<R J J\x-y\>R 

Note that for R > sufficiently large, 

u(x)u(y) | log \x — y\ \ dxdy < — — — / / u(x)u(y) \x — y\ dxdy 



Therefore, 



\x-y\>R R J J\x-y\>R 



$(u)dx + {c + e)\ I j u(x)u(y) log \x - y\ dxdy < F + C(M, 5, R) + C(R, M)M 2 (t) 1/2 

\x—y\<6 



I 2 

By the logarithmic Hardy-Littlewood-Sobolev inequality 

a/2 



/m r 
$>{u)dx-{c + e)— j u\ogudx< F + C(M, 5, R) + C(R, M)M 2 (tY 
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Therefore, choosing k > 1, 

f ui ogu (l^L-( c + e )^]dx+ [ <S>(u)dx<F + C(M,5,R) + C(R,M)M 2 (t) 1/2 . 
J{u>k} \ulogu 2d J J u<k 

By Lemma [8l 

f nlog u ( " (c + e)^l rfx < F + C(M, 5, fl) + C(12, M)M 2 (i). 
J{u>k} \ulogu 2d J 

By LemmaEJ M.2{t) < 1 + t for all t < oo. Since M < M c as defined in (fl4l) . it is possible to choose 
k large enough and e small enough such that $s u>k \ u log udx is bounded on any finite time interval. 
This is sufficient to imply equi-integrability on any finite time interval and hence by Theorem [3] the 
solution u{t) must be global. 

We now refine the argument under the additional hypotheses that tC is bounded below and 
the diffusion is degenerate. First note (see e.g. [5]) that A'(z)z~ 1 dz < oo implies a uniform in 
time bound on the entropy: j &{u)dx > —C(M). Choosing k > 1 this implies with the energy 
dissipation inequality (fTUj) . 



f <$>{u)dx- \ f u1C*udx<FQ + C(M). 

Ju>k ^ J 



Next using the asymptotic expansion of K, at the origin we may choose e and 5 such that 
/ Q(u)dx H ^— / / u{x)u{y) log \x — y\ dxdy 

Ju>k 2 J J {ii>k}n{\x-y\<&} 

(x)\ dx \ kM. 

\J\x\<8 

Finally applying the logarithmic HLS © we have, 



<F + C(M)+ sup \fC(\x\)\M 2 + ( [ |/C(j 
|x|>5 \J|x|<(5 



f ulogu(-^--(c + e)^]dx<F + C(M)+ sup \K{\x\)\M 2 + [ [ \K{x)\dx\kM. 
J{u>k} \ulogu 2d J \ x \ >s \J\x\<5 J 

Therefore by choosing k sufficiently large and e sufficiently small (since M < M c ) we may bound 
j u>k u log udx uniformly in time. By Theorem the solution is global and uniformly bounded. 

We now turn to the proof of (ii), which shows that with additional homogeneity assumptions 
the energy dissipation inequality can be used to deduce optimal decay results for critical problems 
when M < M c . Recall this result is already known |14l lllj . The key is that the scaling invariance 
can be used to transform the energy dissipation inequality into something significantly stronger. 
We proceed by considering the self-similar variables, as in [141 ttH fl3| 13]. defining 9(t, rj) such that 

e- dT 9(T,ri) = u(t,x), (33) 

with coordinates e T r\ = x and e dr — 1 = dt. In these coordinates, if u(t,x) solves ([1]), by the 
homogeneity of the Newtonian potential we have, 

d T 9 = V • (rfi) + A6 2 ~ 2 / d - V • (6VM * 6). (34) 
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Note by definition, uniform boundedness of 6(t) is equivalent the decay stated in (ii) of Theorem 
[5l The key here is that the assumed homogeneity implies that an analogue of the energy dissipation 
inequality (I10p still holds (see for more details). In d > 3, 

— L- J e 2 ~ 2 / d d V + l - J 9(T,r,)\r,\ 2 dr)-± J J 6(r, V )9(t, OM(f] - C)d V d( 

i r Q 2~2/d , , i r a , ,2 , i 



< y c > + 2 y i 7 ?! dr > - 2 y y ^Mcwiv - o^c == s < 

If d > 3 then the argument using the Hardy-Littlewood-Sobolev inequality in [11] applies and 
easily proves the uniform equi-integrability of 9{r,rj). This, in turn, can be shown to imply the 
uniform L°° bound by an easy variant of Theorem [3] (see e.g. [3]). In M 2 , the modified free energy 
becomes, 

Q{6) := j e\og9d n + ^ j \ v \ 2 ed9-^ j j 6^)0(0 tog to -C| c«, 

and the modified energy dissipation inequaliy Q(9{t)) < Q(9q) := Go holds. The key is to use the 
additional second moment term in the modified free energy and Lemma[8]to control the possiblity of 
the entropy or potential energies from being unbounded from below, a trick which was not evidently 
possible in the arguments of Part (i). Using the logarithmic Hardy-Littlewood-Sobolev inequality 

©, 

!"^)/ 0(r,r ] )log9(T, V )dr ] + ^M2(e(T))<go + C(M). 
Assume further of course that M < M c = 8tt. Then by Lemma for all e > 0, 



(9(r,r,)log9(T,n)) + dr, + ^M 2 (e(T)) 



Hence, if we choose e < 2 1 (l — |^) 1 we have, 



(i - ^) / V<r, v ) log 6(t, V )) + d v <g -(i-^j m log ( e ^f j + c(M), 

for all r. Hence 9(t) is uniformly equi-integrable and the result again follows from variants of 
standard continuation arguments. □ 



5 Appendix: Elliptic Estimates 

In this section we discuss basic estimates regarding the elliptic system 

- V • (o(i)Vc) + 7(x)c = /, (35) 

defined on all of W 1 . We assume such estimates, or better ones, can be found elsewhere in the 
classical literature or in textbooks but we could not locate them. The ensuing proofs are technically 
bootstrap arguments, since they require certain norms to be a priori finite to begin with. However, 
the strong solution which vanishes at infinity that we are interested in is constructed by solving (|35p 
on successively larger balls of radius R with zero Dirichlet conditions. On each ball the estimates 
of Lemmas [TUldl] and [12] can be made independent of R and hence can be passed to the limit. We 
skip these straightforward details and only give the formal estimates. 
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Lemma 10. Let c be the strong solution that vanishes at infinity in M. d of 

-V- (a(i)Vc) + 1 (x)c = f, 
for a(x) € L°° strictly positive and^f(x) G L°° non-negative. Ifd>3, then for | + | = |, 1 < q < 
^2 < P < oo we have 



C{p,d) 



(36) 



Moreover, as p — >■ oo, C(p, d) < p. If d> 2, we have the H 1 stability estimate: if f = V • F i/ien, 



l|Vc|| 2 < 



1 



IFIIo. 



(37) 



inf a 

If d > 2 and 7 is a/so strictly positive then for all 1 < p < 00, we have the estimate 

1**5^11/11,. 

Proof. Note this estimate follows from the L p, °° Young's inequality and the representation of c as a 
convolution in the cases when both a and 7 are constants. Define a = (d — 2)p/d — 1 then multiply 
by c a and integrate: 

-/ c «V.„ W V* + / 7W c- (b; = /c«/dx. 
Use inverse chain rule and Sobolev embedding on the LHS (dropping the low order term): 

a J a(x)c a - l \Vc\ 2 dx< J c a fdx 



4a 



(a + 1) 2 



V(C~3~) 



1 



inf a 

11 — 1 1 2 <r 1 

|| C 2 || 2d ^ T-£ — 

5^2 mi a 

| c ||«+i < _L 
1 ~ inf a 



c a fdx 



ia+1 < 



1 



'^("+1) ~ infa" C| W a + 1 ) lu 



dp 



By definition of a, this is in fact 
The second estimate is trivial: 



(inf a)||Vc||| < ||F|| 2 ||Vc||. 

The third estimate is similarly trivial: 

( p-l)/^|Vc K ^ + /7(^-/^ 



(p-l)(^J M* 



Vc^ 2 



dx + (inf7)||c||P < ||c Mp 
1 



p-i 



,C|lp " inf 7 lu 



□ 
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Lemma 11 (Homogeneous gradient estimates). Let d > 3 and let c be a strong solution which 
vanishes at infinity in R rf of 

-V • (a(x)Vc) = /. 

Suppose a(x) G C 1 is strictly positive, bounded and Va is a/so uniformly bounded. Then if p > 
dj (d — 2) we have, 

\\D 2 c\\ p <C(a,p,d)(\\f\\ p + 



2p+d 

Moreover, as p — > oo, C(a,p,d) < p. 

Remark 10. It is crucial to note that 1 < pd/ (2p + d) < p. 

Proof. The proof is a variant of Theorem 9.11 in [23] adapted to IR^. The following Gagliardo- 
Nirenberg-type inequality, (see e.g. Theorem 7.28 in [23]), will be used: 

l|Vff|| P < d , P hlff 2 \\D 2 g\\l /2 < 4D 2 g\\ P + -J\g\\ P . (38) 

Consider a ball of radius R > centered at a point xo £ denoted Br(xq). By the Calderon- 
Zygmund inequality we have for a sufficiently smooth v supported in Br(xq) and 1 < p < oo, 

II^HlpCBb^o)) - °( a; o)||At;||p. 

Hence, 

II^IIp £ l|V • ((o(x ) - a(x))V«)||p + ||V • (o(s)Vt;)||p. 
We control the first term using l38| for some e > 0, 

||V • (a(xo) — a(x))Vv\\ p < sup |a(xo) — a(x)| \\ Av\\ p + || Va • S7v\\ p 

x&B R (xo) 

< sup |a(xo) — a(x)\ \\D 2 v \\ p + ||Va||oo||Vu || p 

x€B R (xo) 



< sup |a(xo) — a(x)| ||Z? 2 ti|| p + HValloo ( e||D 2 w|| p H 

x£B R (x ) V e J 



x£B R (x ) 

Hence, choosing R sufficiently small, depending on the smoothness of a(x) and choosing e sufficiently 
small, we can then deduce, 

\\D 2 v\\ p <\\V-(a(x)Vv)\\ p + \\v\\ p . 
Choosing v = r\c in the above inequality gives, for any a £ (0, 1), 

\\D 2 c\\ Lp{B{tr{xo)) < ||V • a{x)Vc\\ LP{BrTRixo)) + \\aVrj ■ Vc 

+ caArj + cVa • Vrj\\ p + ||Vo|| L cx, (B(7h) ||c|| L p( B(7K ) 

c c c 

< \\f\\LP(B aR (x )) + ^ll^VcHip^^)) + -jp\\ca\\ L p { B aR (xo)) + ^ll cVa ll^(B CT fl(x )) 
+ WV a \\L°°(B< rR {xo))\\ c \\LP(B< rR (xo))- 

We now have to deal with the latter error terms. Using the interpolation inequality on page 250 of 
[23], which is essentially just (|38h . we have, 



C 

ll«Vc|| L p (B(Tfl ( x . o) ) < ||a||oo||Vc|| L p (B(Tfl(:!;o )) < ||a|| TO ife||.D 2 c|| 1 ^( Baii ( X0 )) + \\a\\oo^\\c\\LP(B aR (xo))> 
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where crucially the constant C does not depend on a or R (C is some power of a constant of a 
Gagliardo-Nirenberg inequality which depends only on the geometry and not the diameter). Hence 
by choosing e sufficiently small independent of R we have, 

\\D 2 c\\ Lp{Bitr{xq)) < \\f\\ L p(B aR (x )) + (Jp + ^ + \\Va\\ L °o {BaR{xo)) J \\c\\ LP ( BaR {x ))- (39) 

Since p > d/(d — 2), (j36[) ensures that latter term will add up as we cover M. d with balls (using that 
we may take R > 1 due to the bound on Va). Therefore, 

\\D 2 c\\ p <\\f\\ p + \\f\\^. 

2p+d 

□ 

Lemma 12. Let d > 2 and in ]R d let c be a strong solution which vanishes at infinity in W 1 of 

-V- (a(i)Vc) + 7 (x)c = /, 

for-)(x) G L°° strictly positive and a(x) G C 1 strictly positve, bounded andVa is uniformly bounded. 
Then we have the gradient estimate, 

\\D 2 c\\ p <C(a,p,d)\\f\\ p , 

with C(a,p,d) < p as p — )• oo . 

Proof Easy variant of the above using the better L p norm estimate (|37p to control the lower order 
terms in (1391), □ 



6 Appendix: Gagliardo-Nirenberg 

Lemma 13 (Homogeneous Gagliardo-Nirenberg). Let d > 2 and f : R d -)• R saiis/y f £ L p n L q 
and V/ fc € L r . Moreover let 1 < p < rk < dk, k < q < rkd/(d — r) and 

Iks ... 

< 0. 40 

r q d 

Then there exists a constant Cqns which depends on s,p,q,r,d such that 

||/lk* < C GN s\\f\\ a L l\\f k \\^ s , r , (41) 

where < oti satisfy 

l = ai k + a 2 , (42) 

and 

= ai(— + • 43 

q p d r p 
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